Out of plane screening and dipolar interactions in heterostructures 
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Out-of-plane screening (OPS) is expected to occur generally in metal-semiconductor interfaces 
but this aspect has been overlooked in previous studies. In this paper we study the effect of 
OPS in electron-hole bilayer (EHBL) systems. The validity of the dipolar interaction induced 
by OPS is justified with a RPA calculation. Effect of OPS in electron-hole liquid with close-by 
screening layers is studied. We find that OPS affects the electronic properties in low density and 
long wavelength regime. The corresponding zero-temperature phase diagram is obtained within a 
mean field treatment. We argue that our result is in general relevant to other heterostrucutures. 
The case of strongly correlated EHBL is also discussed. 

PACS numbers: 77.80.bn, 71.35.Ee, 71.35.-y 



I. INTRODUCTION 

Modern micro-electronics relies to a large degree on 
surface science, which concerns the material properties 
near a surface or interface. To enhance the performance 
of such devices, knowledge of the electronic states near 
the interfaces is required. Near a surface or interface, 
electronic reconstruction may alter three key factors - 
interaction strengths, bandwidths and electron densities 
[l[ which determine electronic states and their properties. 

In this paper, we consider another factor - the mod- 
ification in form of interaction between electrons. For 
instance, in an insulator-semiconductor-insulator super- 
structure, if the dielectric constant of the semiconductor 
is sizably larger than that of insulator (barrier layer), the 
image charges induced at semiconductor-insulator inter- 
face can substantially enhance the binding energy of the 
excitons confined in the semiconductor layer @, I n 
this case, the electrons and holes do not interact via usual 
Coulomb potential after the effect of the image charges 
at the semiconductor-insulator interface is taken into ac- 
count. 

Recently, Huang et al. observed non-activated elec- 
tronic conductivity of a two-dimensional (2D) low den- 
sity hole system in a hetero junction insulated-gate field- 
effect transistor Such non-activated conductivity is 
unexpected as at low charge density strong Coulomb in- 
teraction is expected to crystallize the system (Wigner 
crystal), which is then pinned by disorder resulting in in- 
sulating behavior and activated conductivity. Huang et 
al. attribute the behavior to the screening of Coulomb 
interactions by the metallic gate, which leads to destruc- 
tion of the Wigner crystal phase. Physically, the metallic 
gate which is located at a distance away from the 2D hole 
gas, provides an out of plane screening (OPS) to the hole- 
hole interaction, resulting in effective dipolar interaction 
between holes. Microscopically, when a charge is placed 
near a metal surface, an image charge of opposite sign will 
be induced at the surface to screen out the (static) elec- 
tric field from the charge. From elementary electrostat- 
ics, the system can be described equivalently as a dipole 
formed by the charge and its image charge and the in- 



teraction between two charges located near the interface 
changes from a Coulomb potential ~ 1/r to a dipolar 
potential ~ 1/r 3 . This modified interaction, which is 
generally expected to exist in metal-semiconductor het- 
erostructures, can change the electronic properties near 
the interface. Unexpectedly, there has been no detailed 
theoretical study of this effect on electronic properties 
until recently [a]. The neglect of OPS might be due to 
dynamical screening of in-plane charges J5|. For high 
charge density, the screening can effectively reduce both 
Coulomb and dipolar interactions to short range interac- 
tions. However for low charge density electronic liquids 
in-plane screening is less effective and OPS can lead to a 
difference, as is observed by Huang et al. Q. 

In this paper, we study how OPS affects the electronic 
properties in systems with two-layer of charges of oppo- 
site sign, i.e. the 2D electron-hole bilayer (EHBL) sys- 
tem. We shall study how OPS affects Wigner crystaliza- 
tion and exciton condensate in the system 0, @ and will 
also comment on the effect of OPS in interfaces between 
metals and strongly correlated electron systems 0-0- 



II. OPS AND EFFECTIVE INTERACTION 
BETWEEN CHARGES 

In this section we provide the details for the EHBL 
systems we study and the corresponding OPS effective 
interaction. We shall assume that the only effect of the 
metallic screening layers is to provide an image charge 
for point charges sitting close to it and the effective in- 
teraction between charges will be derived from the im- 
age charge picture. The validity of this approximation is 

(s) 

bounded by the plasma frequency ujp of the screening 
layer, above which the screening layer cannot respond 
rapidly to the charge fluctuations. Thus our approxima- 
tion is valid when the plasma frequency of the EHBL 

(s) 

layer tu p is much less than uip , or that the screening 
layer has density of electric charge much larger than the 
charge density of the EHBL layers we consider. The im- 
age charge picture can be justified by a Random Phase 
Approximation (RPA) calculation which is shown in the 
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FIG. 1: (a) EHBL system separated by distance b. (b) EHBL 
with OPS by metallic plates in both layers. Dotted line rep- 
resents metallic interface, separated from the main layer by 
a distance of a/2, (c) Similar to (b) but with only one OPS 
layer, (d) Charges (black dots) and screening charge response 
at the metal interface (grey patches), (e) Effective image 
charge (grey dots) and effective interactions 1 [/ mtra and V x . 
(f) Charges attract when they are aligned while repel when 
they are not. This behavior is different from the Coulomb 
potential which is always attractive for a pair of electron and 
hole. The repulsive behavior inhibits exciton pairing. 



Appendix. 

Starting with a EHBL system (Fig.QJa)), two metallic 
screening layers can be added as shown in Fig.QJb), or a 
single metallic screening layer can be added as shown 
in Fig. [lie). We first consider the two-layer case (b). 
Fig. [lid) depicts the charge response in the metallic layer 
to a nearby charge. The charge response is assumed to 
be an image charge, which carries opposite charge of the 
same magnitude and is centered at distance a from the 
point charge. Thus the point charge and the screening 
charge together form a dipole. We have assumed that 
the distance between the two layers of charges b is suffi- 
ciently larger than a (b ^> a) such that the presence of 
the other screening layer does not affect the simple dipole 
picture. In this case, the intralayer interaction between 
two charges located in an OPS layer (Fig. (He)) is in real 
space 
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where r is the charge-charge distance within the charge 
plane. 

It is easy to see that for r ^> a, y mtra scales as 1/r 3 
while for r <C a it follows the usual Coulomb scaling 1/r. 

By using 2D Fourier transform 
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For an electron and a hole sitting in different layers, the 
interlayer interaction is 
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and its Fourier counterpart is 
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t Et h and e x are the intra-layer and inter-layer dielectric 
constants, respectively. 

Next we consider EHBL with only one metallic screen- 
ing layer (see Fig. [He)). In this case the two layers of 
charges have distance a/2 + 6 (layer 1) and a/2 (layer 
2) from the screening layer, respectively. The intralayer 
interactions are thus 



e 2 1 

V; nt ™(r) = — - 



1 



V 2 a (r) = - - 



ei \r yjr 2 + (a + 2b) 2 
2 n 1 



02 



\/r 2 +. 



with corresponding Fourier transforms 
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The corresponding intralayer interaction is given by 
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III. COLLECTIVE DENSITY RESPONSES 

In this section we study the collective density responses 
of the EHBL systems we considered. For a two compo- 
nent electronic system, the density-density response of 
the system is described by a 2 x 2 matrix Xij(q,w) with 
i,j = 1,2. The density-density response matrix is given 
in RPA by 
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where 
n{q,u) 



(6) 



(1 - X0l(<7,w)Vii(q))(l - X02(g,w)V22(g)) 

-xoi(g,w)Pi2(g)xo2(g,w)V2i(g) , (7) 



Vij- (g) is the "bare" interaction between i and j com- 
ponents of the electronic liquid and 

d 2 k H F ( £ i 8) ) - n F (4+J 



xoi(g,w) 



where e£ ~ k 2 /2m^ is kinetic energy of species i par- 
ticles (of mass m'*'), iif is the Fermi-Dirac distribution 
function and g s — 2 is spin degeneracy. In the case of two 
screening layers the interactions Vii(22) an d V\2 = V21 are 
given by l/ intra (q) (eq. ©) and ^(g) (eq. ©), respec- 
tively whereas they are given by V"^ a (q) (eq. (g))) and 

F s,1 (g) (eq. (JSJ), respectively if there is only one screen- 
ing layer. 

Next we study the collective excitations (i.e. plasmons) 
in the system. The dispersion of the collective excitations 
are given by the equation 
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K(q,uj{q)) = . 



(9) 



We shall first consider the long wavelength limit (q —> 0) 
where the equation can be studied analytically. In this 
limit it is easy to show that 
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where n 



is carrier density. We have neglected 
the component index i for brevity. 

We begin with the Coulomb case (no screening layer). 

The interactions are respectively Vi t 2{q) = 7^ and 



equation in q - 
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We first consider the case e 2 = t\e 2 such that the term 
in the square bracket is zero. In this case we need to 
expand the interlayer interaction to one order higher in 

q. As a result the last term in eq. (JTTJ) is replaced by 

3 

a term of order £L and the plasmon equation at long 
wavelength limit yields two solutions, which are the out- 
of-phase mode (ui ~ q) and in-phase mode (lu ~ */g). In- 
deed this occurs usually in a 2D electronic systems with 
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FIG. 2: Plasmon excitations in EHBL with two screening 
layers computed numerically. The spectrum is computed by 
solving eq. ([9]) numerically with the full expression of xo(9,w) 
flol |. The solid lines are the plasmons excitations and the 
dash line represents the boundary of the particle-hole contin- 
uum. The wave-number q and frequency w are normalized 
with respect to Fermi momentum &f and Fermi energy tF, 
respectively. We set a = 1, 6 = 15, mi = 1, m2 = 1.25, 
kF = 10 and all e = 1 in the calculation. 



both conduction and valence bands where the same di- 
electric constant e 2 = t\t 2 is found for all interactions. In 
the more general case e 2 ^ £162, which arises quite natu- 
rally in the complex environment of EHBL hctcrostruc- 
tures, we can easily see from eq. (fTTj) that the plasmon 
frequency scales as w ~ ^fq. There are two modes of 
plasmons. 

For the OPS case with two screening layers, the in- 
teractions are respectively V\^ 2 = fff^(l ' 
2^e 2 / ei , 2 (a - a 2 g/2) and V x = -^V« b (l 
— 2na 2 e 2 q/ e x for q <C a -1 . Notice the removal of the 
1/q singularity in the interactions by OPS. We then ob- 
tain after solving the equation the collective modes (up 
to order q 2 ) 
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Notice that OPS effectively reduced the long-ranged 
Coulomb interaction into short-ranged interactions re- 
sulting in two collective modes scaling linearly with q. 
The collective modes represent separate collective mo- 
tion of the two layers because (V x ) 2 is of higher order 
in q than V\V 2 , and the inter-layer interaction appears 
only to order q 3 . For completeness, we have computed 
numerically the collective modes spectrums at finite q as 
shown in Fig. [51 

With only one screening layer, the interactions 
are Vi(q) ~ ((a + 2b) - (a + 2b) 2 q/2) , V 2 {q) ~ 

^(a-a 2 g/2) and V x (q) = -2l^L e ~ib(i _ e -ga) _ 
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FIG. 3: Plasmon excitations in EHBL with only one screening 
layer. At small q, both plasmon modes scale linearly with q 
with a larger slope oc v^6 + a for u>%. We set a = 1, b = 15, 
mi = 1, m,2 = 1.25, kF = 10 and all e = 1 in the calculation. 



— 2Trae 2 /e x , respectively at small q. The collective modes 
are given by (up to order q 2 ) 
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Again there are two linear plasmon modes and effect of 
V x does not enter until q 3 . The main difference is that 
the electron-hole layer separation b enters the slope of u>i 
mode (oc \j2b + a). The numerically calculated plasmon 
spcctrums are depicted in Fig.[3l 



IV. EXCITON CONDENSATION AND WIGNER 
CRYSTALIZATION 



In this section we study exciton condensation and 
Wigner crystalization in an electron-hole liquid with 
OPS. The system without OPS has been extensively 
studied for the search of exciton condensation. We shall 
consider exciton condensation in a BCS type mean-field 
theory where the exciton condensation is described by 
the order parameter (cikfC 2k ^) (1,2 are layer indices). 
For simplicity we assume the layers are doped with equal 
amount of charges (with opposite signs) and the electrons 
and holes are spin-polarized. Singlet pairing of excitons 
is implicitly assumed. 

The EHBL Hamiltonian in momentum representation 
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where a = 1,2 is the layer index; c k (ct) is the momen- 
tum k fermion annihilation (creation) operator, = 

2^ Ha is the electron or hole dispersion and V a (k) 

(V x (k)) is the intralayer (interlayer) OPS effective in- 
teraction. Next we employ the standard Hartree-Fock- 
Bogoliubov method [ll[ to derive the mean field equa- 
tions for exciton condensate. The Hartree-Fock terms 



s fc = J2 q va (P - <?) \c ak c a kj modify the particle dis- 
persions £j? — > — and need to be solved self- 
consistently. Here we concentrate on the effect of ex- 
citon binding on the Fermi surface and shall assume 
that the self-energy can be captured by introducing effec- 
tive masses to* (e a ) and renormalized chemical potentials 
/i* (e Q ), i.e. f£ - Eg ~ g^F ~ Ma- With this approxima- 
tion, we obtain the mean field Bogoliubov Hamiltonian 
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is the exciton order parameter. Hmf can be diagonalized 
easily by the Bogoliubov transformation 
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(to*" 1 + m* 2 ~ L )/2 and fj, = + f4)/2. 



The ground state wavefunction is 

\i>a) = I] Uk + Vkc\ k c\A |0) 



(19) 

where Afc is determined by the self-consistent equation 

^ = -\T, v ^ k -i)^-- ( 2 °) 

9 q 

The equation is to be solved with the particle number 
constraint 
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FIG. 4: The phase diagram for EHBL with two- and one- 
screening layers for varied image charge separation a. b is the 
bilayer separation and r s ~ is the dimensionless average 
particle-particle separation in unit of the effective Bohr radius 
as = € x h 2 /m e ffe 2 . Below the transition lines b c (r s ) an exciton 
gap A of magnitude larger than 10~ fj, is formed. 



where v\ is the probability of finding an electron-hole 
pair in state k at the ground state. A zero-temperature 
phase diagram can be determined by numerically solving 
eqs. (gni) and (JH]). 

To simplify calculation we assume further that exciton 
gap is momentum independent A^ = A and A is deter- 
mined by minimizing the ground state energy. We note 
that we are considering a band structure with isotropic 
dispersion and the electron and hole Fermi surfaces are 
perfectly nested. In this case, the exciton pairing gap 
A is always non-zero in the mean-field theory, although 
its value can be very small. In reality the mean-field gap 
will be destroyed by quantum fluctuations when it's mag- 
nitude is small [12], but this is not reflected in a mean 
field theory. To capture this physics qualitatively, we 
assume that the transition from the exciton condensed 
state to the normal state occurs at A = 10~ 5 /i. Al- 
though quantitatively unreliable, this procedure allows 
us to examine the effect of screening on the phase dia- 
gram semi-quantitatively as we shall see below. 

With the above criteria, the phase diagram for dif- 
ferent average particle-particle separation r s = -^\J~-^[ 

(n is particle/hole density; as — e x h 2 /m c ffe 2 is the ef- 
fective Bohr radius of electron-hole pair) and transition 
layer separation b c (r s ) can be determined by solving the 
self-consistent equations (|20|) and ([2l"j) . We first consider 
EHBL with two screening layers. The result of calcula- 
tion is depicted in Fig. |4] for different separation between 
the electron/hole and its image charge a (filled symbols). 

In the small r s (high density) regime, kinetic energy 
dominates over potential energy and the exciton pairing 
gap goes to zero as r s — > 0. In large r s or low density 



limit, the exciton pairing is diminished due to the repul- 
sive nature of the interlayer OPS potential at short dis- 
tance (see Fig.QJf)). This leads to a linear dependence of 
V x (k) versus k at small k (see eq. (O). In this case, the 
gap equation eq. (|2T)]) is of the form f Q F —j===d 2 k 

cons' an: for small gap A, where kp ~ l/r s and larger 
r s (smaller kp) implies a smaller A to satisfy the equa- 
tion. The electrons and holes need to be placed closer to 
each other to produce a large enough A and leads to the 
drop of b c (r s ) at large r s . As the screening separation a 
increases, the transition line shifts upward as the inter- 
layer OPS potential is strengthened which enhances pair- 
ing. For a = 25 (comparable with b), the image charge 
effect becomes negligible and potential becomes essen- 
tially Coulomb-like which permits exciton formation for 
all r s we considered (cf. Fig. 1 in Ref. fl2j). The main 
effect of OPS potential is to suppress exciton pairing at 
low density. 

Previous numerical study of the same EHBL with 
no screening layer [l2j reveals also an excitonic Wigner 
crystal phase at large r s . Wigner crystal is commonly 
formed in low density (i.e. large r s ) electron liquid be- 
cause of domination of Coulomb repulsive potential en- 
ergy (~ l/j's) over kinetic energy (~ 1/r 2 ). To minimize 
the potential energy the electron wavefunction "crystal- 
lizes" to ensure maximum separation between electrons 
which yields the Wigner crystal phase. Here we argue 
that OPS suppresses the Wigner crystal phase in two 
ways. Firstly, as shown above, exciton formation is sup- 
pressed at large r s and thus the excitonic Wigner crys- 
tal is unlikely to form. On the other hand, electronic 
Wigner crystals in separated layers are also prohibited 
since introduction of OPS reduces the (intralayer) poten- 
tial energy and changes its scaling form to ~ 1 /r|? (dipo- 
lar interaction, see eq. (|A.1[1 ) at large particle separation 
r 3> a. In this case kinetic energy again dominates at 
large r s and an usual electron/hole liquid phase should 
occur. The situation is similar to the case as found in 
Ref. [J] where the electronic Wigner crystal phase is de- 
stroyed by screening. We note, however that our simple 
study cannot rule out the possibility of having a Wigner 
crystal phase at some intermediate values of r s where the 
kinetic and potential energies are of comparable magni- 
tudes. 

We now consider the situation of EHBL with only one 
screening layer which may be easier to realize experi- 
mentally (Fig.OJc)). In this case we adopt eq. J5]) for 
interlayer interaction, where V x (k) scales as constant at 
small k. We can again consider the gap equation and ar- 
gue similarly that the exciton phase boundary would also 
drop at large r s , as in the two- layer screening case. In- 
deed we have solved the gap equations and find that the 
phase diagram is qualitatively the same as the two OPS 
layer case except that the area under the phase boundary 
b c (r s ) is larger (see Fig. @] (open symbols)). 

For the Wigner crystal phase, the "asymmetric" OPS 
introduces some complications. First we note that an 
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excitonic Wigner crystal phase is also unlikely to occur 
at large r s . However the system may form a hybrid 
phase where a Wigner crystal is formed at layer 1 and 
electron/hole liquid phase remains for layer 2 because 
screening mainly affects layer 2. To examine this pos- 
sibility we check the effective intralayer interaction af- 
ter taking into account the screening effect of the other 
charged layer (see eq. (|A.1[) in Appendix and discussions 
thereafter). We see that the effective intralayer inter- 
action is mainly dominated by V / 1 ™ tra (g), and screening 
from the other layer is not important. Therefore, we 
expect that at large r s kinetic energy again dominates 
and the both layers are in the electron/hole liquid phase. 
Notice, however that V r 1 lntra (g) has a dipolar form only 
when r s ~ r/as 3> b/as for layer 1. Thus for some large 
enough h/as, a hybrid phase (Wigner crystal at layer 1, 
electron/hole liquid at layer 2) may still occur at some 
intermediate densities &/ob > r s > 1. 

We see that OPS becomes important for low density 
electronic systems due to change in scaling of the poten- 
tial energy. Generally speaking, for heterostructures, in- 
sulating behavior resulting from low carrier density can 
be avoided by addition of metallic screening layers Q. 
This method may be preferred over other methods like 
increasing carrier density by dopants since dopants act 
like impurities and introduce unnecessary scattering at 
low temperature. 



V. STRONGLY CORRELATED EHBL 

In strongly correlated materials, the basic electronic 
properties are determined by the bandwidth, the on-site 
Coulomb interactions U and the charge transfer energy 
E c . If such a ultra-thin film, originally a Mott insula- 
tor, is placed close to a metal surface, U and E c can be 
strongly reduced by OPS Q. When the bandwidth ex- 
ceeds the suppressed U and E c , the insulating film can 
undergo an insulator-metal phase transition. Further- 
more, if a heterostructure is formed, structural relaxation 
and local electronic states may exist at the interfaces. For 
instance, in an interface formed by YBa2Cus07 (YBCO) 
cuprate and metal [8, 9], the Cu02 plane near the inter- 
face (depletion layer) is intrinsically doped by electronic 
reconstruction resulting in a strongly correlated electron 
system with OPS interaction induced by the metal. We 
shall consider here how OPS would affect the properties 
of this system. 

The mean field analysis on effect of OPS can also be 
performed for strongly correlated EHBL systems [HI [l4| 
with a two-layer t-J type model. We assume here that 
the suppression of U and E c induced by OPS are not 
strong enough to destroy strong correlation, otherwise 
we can simply apply the usual electron-hole liquid pic- 
ture described in previous section. Therefore the set- 
ting is similar to that shown in Fig.QJb) except that 
the electron-hole liquid is replaced by a strongly corre- 
lated EHBL with holons and doublons and the excitons 



are formed by holon-doublon pairs instead of electron- 
hole pairs. A mean field calculation similar to that 
of Ref. [HI can be carried out by applying the slave- 
boson mean field theory to the two-layer t-J model. The 
main difference is that the on-site interlayer interaction 
J2i blibubl^i is replaced by the OPS effective inter- 
action J2ij Vfiblibublfoj, where b ai {tf ai ) is the bosonic 
holon (a — I) or doublon (a — 2) annihilation (creation) 
operator of layer a at site i. The OPS interaction is then 
decoupled as 



ij pqk 

« J2K(Kb2 P + bt p b\ p ) 

P 



where A^ = Y^ q Vp-q (bi q b2q) is the exciton pairing. As- 
suming that Ap = A b S(p) is homogeneous in space, we 
obtain a mean field Hamiltonian which is of the same 
form as in previous study [HI for on-site interaction Vb 
with (&ii&2i) ~ J2k (kife&2fc)- Since in terms of exciton 
pairing the attract-repel behavior renders the interlayer 
OPS interaction resembling an on-site interaction (see 
Fig.QJf)) , we expect that the mean field phase diagrams 
in both case are qualitatively the same. The introduction 
of OPS interaction solely shifts the exciton phase bound- 
ary due to a reduction of interaction strength, as in the 
case of usual electron-hole liquid. 

We next comment on the possibility of forming spa- 
tially inhomogeneous phases. One example of such inho- 
mogeneity is charge corrugation in the form of stripes. 
By applying mean field theory to t-J model with long 
range Coulomb interaction VcY^tyj -f-. n i n ji it is shown 
that stripes arc preferred to minimize the exchange J 
term [15| . In particular, it is the decoupling of the ex- 
change term into the anti-ferromagnetic channel that 
drives the stripe formation, while the Coulomb interac- 
tion controls the spacing evolution of stripes with dop- 
ing. Moreover, the stripes spacing increases as the dop- 
ing S decreases. The effect of OPS on stripes is two- 
fold. Firstly, OPS weakens the on-site repulsion U [l6j] 
and thus the superexchange J ~ t 2 /U term is enhanced 
(assuming that strong correlation is still intact). Con- 
sequently, the stripes phase is strengthened. On the 
other hand, Coulomb interaction tends to smooth out 
the charge density, while a dipolar interaction (V ~ 1/r 3 
for large r) would be less effective and a more inhomo- 
geneous phase would be preferred. Notice that extreme 
charge inhomogeneity like phase separation 17] is not 
likely since the OPS interaction scales like 1/r for small 
r and still suppresses phase separation. 
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VI. SUMMARY 

We have constructed a dipolar interaction for OPS ef- 
fect of metallic layer in heterostructures and have justi- 
fied the construction by a RPA calculation. The OPS 
interaction is expected to be present rather generally at 
interfaces with metallic layers. We apply the OPS inter- 
action to EHBL system and find that OPS mainly affects 
the electronic properties in the low density regime. Our 
conclusion is not restricted to EHBL since the behavior is 
mainly due to the modification of the interaction scaling 
from 1/r (Coulomb) to 1/r 3 (dipole) at distance of large 
r. OPS might be employed to eliminate Wigner- crystal 
like behavior at low temperatures. For strongly corre- 
lated electron systems, OPS mainly affects the magnetic 
channel by reducing the Hubbard U and charge transfer 
energy E c . The reduction of U may drive the system 
into usual electron liquid. Furthermore, the reduction in 
interaction range may drive the system into an inhomo- 
geneous state. 



2_2 
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XOs 
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FIG. 5: Diagram for construction of OPS interactions. Thin 
lines are bare interactions, bubble is the density-density re- 
sponse function XOs and the thick line is the resulting effective 
interaction. 



therefore 

2Tre 2 N F 
F ° = q + 2ve 2 N F ~ 1 ' 

This approximation is valid if the charge density n 2 of 
layer 2 is much less than the density of the screening 
layer n s and q < yfni <C We shall take the same 

limit in the following derivations. This gives eq. @. 

Similarly we can construct the interlayer interaction 
eq. ©: 
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Appendix: Justification of OPS interactions by RPA 

In this appendix, we employ RPA to justify the image- 
charge picture of OPS interactions. The RPA method 
enables us to obtain an effective interaction by "integrat- 
ing out" the screening layers. 

First we consider a charged layer 2 with a metallic 
screening layer s separated from layer 2 by distance a/2, 
as shown in Fig.QJc). We can write down the effective 
intralayer interaction of layer 2 after taking into account 
the effect of screening by the metallic layer (see Fig. [5]) 



V mtIa (q) = V 2 {q) + 



V2s(q)xosVs2(q) 
1 - XosV S s(q) 



= 2ne z - - 



2ne 2 



1 e- aq 2Tre 2 N F 



q q + 2ne 2 N F 



(1- 



where V 2 (q) = V ss (q) = 2ite 2 jq are the bare Coulomb 
interactions of layer 2 and screening layer s, V2 S {q) = 
Vs2(q) = 2-Ke 2 e~ aq / 2 /q is the interlayer Coulomb interac- 
tion between the layers, and xos = Xos{q — > 0,w = 0) = 
— N F is the (7 — ^ static density-density response func- 
tion of layer s, N F is density of states at the Fermi 
surface. Notice we have assumed that q is small (long 
wavelength limit) in writing down the interactions and 



Vis{q)xQsVs2{q) 

1 - X0sV ss (q) 



2ne 2 



,-bq e ~(a+b)q 



-Fr 



q 



p-bq 
27TC 2 — (1 



q 

-aq\ 



where V 12 (q) = 2Tre 2 e- b y q , V ls (q) = 2ne 2 e- ( - a / 2+b ^ /q 
and V S 2(q) = 2ne 2 e~ aq / 2 / q are the bare interlayer inter- 
actions between the pair of layers (1,2), (l,s) and (s, 2) 
respectively. For two screening layers (Fig.QJb)), we as- 
sume that layer 2 and s form an effective system 2' and 
thus we can adopt V x ' 1 as the "bare" interlayer interac- 
tions in the following: 



V^(q) = V*>\q) + 



, , V ls (q)xosV s2 >{q) 



2Tre z l~e- aq 



1 - XosV ss (q) 

e -bq e -{a+b)q 



-Ft 



q 



q 



2ne 2 



where F s2 '0?) = V x > \b 



-bq 



(l-e- aq ) 2 , 



-(i 



+b)q M 



-aq\ 



1 21 g 

The validity of assuming the effective system 2' is based 
on the choice of b ^> a. This gives eq. ([3]). 

Here we derive the effective intralayer interaction 
V[ffi(q) with two-layer OPS (Fig.QJb)) taking into ac- 
count the screening of system 2' (i.e. integrated out all 
the screening by 2, si and s2): 



intra 



(q) 



X02 (V x ' 2 (q)) 2 

1 _ V intr a{q)xQ2 



(A.l) 
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In the small q limit, V x ' 2 and V r 2 lntra scale as q and con- 
stant respectively. The second term due to screening is 
of higher order in q and thus it cannot alter the scaling 
of the V 1 lntra (g) term (~ constant). One can repeat the 



analysis for the one-layer OPS case (see Fig.fljc)) and the 
scaling of the effective intralayer interaction Vf^ tra (g) in 
the lowest order of q is not affected by screening of the 
opposite charged layer. 
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